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Abstract — Coupled elliptic quantum dots with different as- 
pect ratios containing up to two electrons are studied using a 
model confinement potential in the presence of magnetic fields. 
Single and two particle Schrodinger equations are solved using 
numerical exact diagonolization to obtain the exchange energy 
and chemical potentials. As the ratio between the confinement 
strengths in directions perpendicular and parallel to the coupling 
direction of the double dots increases, the exchange energy at 
zero magnetic field increases, while the magnetic field of the 
singlet-triplet transition decreases. By investigating the charge 
stability diagram, we find that as inter-dot detuning increases, 
the absolute value of the exchange energy increases superlinearly 
followed by saturation. This behavior is attributed to the electron 
density differences between the singlet and triplet states in the 
assymetric quantum dot systems. 

Index Terms — Simulation; quantum dots; exchange energy; 
exact diagonalization. 

I. Introduction 

Coupled quantum dots (QDs) are of particular interest for 
spin-based quantum computation because quantum logic gates 
(such as a Controlled-NOT gate) can be realized via the 
interaction between two spin quantum bits (qubits), i.e., the 
spins of two electrons, each trapped in an individual quantum 
dot [1]. In such devices, the interaction between the two spins 
is proportional to the exchange energy J, which is equivalent 
to the splitting between the lowest singlet and triplet two- 
electron states. 

Extensive theoretical and experimental works have been 
done to study the exchange interaction in coupled semiconduc- 
tor QD systems. From a theoretical point of view, variational 
methods such as the Heitler-London method or exact diago- 
nalization are commonly used to obtain the dependence of the 
exchange interaction on the system parameters, e.g., the inter- 
dot separation, the tunneling barrier between the QDs, and the 
externally applied magnetic field [2]-[5], Experimentally, the 
extraction of the exchange interaction relies on the charge sta- 
bility diagram in which the boundaries between distinctive sta- 
ble charge states, i.e., between the states with fixed number of 
electrons (N±, N2) in each of the coupled dots, are represented 
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as functions of the two controlling gate biases, one for each dot 
[6]. To obtain J, two different methods were used: one is based 
on the Hubbard model, which involves the analysis of double- 
triplet point (DTP) separation on the stability diagram [7]; the 
other utilizes Zeeman splitting to measure the exchange energy 
as a function of inter-dot detuning in a coherent control cycle 
of state preparation, spin-interaction, and projective readout 
in laterally coupled semiconductor quantum dot systems [8]. 
In such an experiment, the charge state of the two qubits is 
manipulated through a charging cycle from the (0, 2) state to 
the (1, 1) state and finally back to the (0, 2) state in a coherent 
way. These kinds of experiments are stimulating interest in the 
investigation of the exchange interaction and the associated 
charge re-distribution as a function of inter-dot bias detuning 
in coupled QDs. 

In experimental planar laterally coupled quantum dot device 
structures, the strength of the confining potentials in the 
coupling direction and the direction perpendicular to it in the 
plane of the two-dimensional electron gas are different due 
to the top gate patterning [9], [10]. In this paper, we account 
for this effect and study coupled quantum dots with different 
aspect ratios in the confinement potential, i.e., with different 
ratios of the confinement strengths in the two directions. 

We introduce our model Gaussian-shaped confinement po- 
tential and exact diagonalization method in Section II. In 
Section III, we discuss the different dependences of the 
exchange interaction on the magnetic fields for different dot 
deformations. We then analyze the stability diagram for var- 
ious coupled QD systems followed by a comparison of the 
dependences of the exchange energy on inter-dot detuning 
for dots with different aspect ratios and their relationship to 
electron density localization. Finally, we present concluding 
remarks in Section IV. 

II. Model and Method 

We use the following model potential to describe the cou- 
pled quantum dot system [5]: 

V(r) = -V L e- {x+d/2)2/R2 * +v2/R2 »-V R e- {x ~ d/2)2/R * +v2/R2 y, 

(1) 

where Vl and Vr are the depth of the left and right dots 
(equivalent to the QD gate voltages in experimental structures 
[6]) which can be independently varied, d is the inter-dot 
separation, R x and R y are the extension of the each dot in 
the x and y direction, respectively. 

Our computational approach consists of two steps [11]. In 
the first step, we solve a single-particle problem with the 
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Hamiltonian 



Here, 



h(r) = 



^t(p + -A) 2 

2m c 



V(r) 



(2) 



= 0.067m e is the electron effective mass in 
GaAs, and A = |[— By, Bx] is the vector potential of the 
magnetic field B oriented perpendicular to the xy-plane. This 
Hamiltonian is diagonalized in the basis of the product of 
two harmonic oscillator states in each direction, ipi(r) = 
4>n(x)(f>m.(y), where <p n (<fr m ) denotes the n-th (m-th) har- 
monic oscillator states (in this work, we use eight harmonic 
states in each direction). 

In the second step, we solve the two-particle problem for 
which the Hamiltonian is given by 



where C(ri,r 2 ) 



tf(r 1 ,r 2 ) = /i(r 1 ) + /i(r 2 )+C(r 1 ,r 2 ) 
e 2 1 



(3) 



accounts for the Coulomb 



e |ri - r 2 

interaction between the two electrons and e = 12.9 is GaAs 
dieletric constant. The Zeeman effect is not included in (3) 
since its effect trivially lowers the energy of the triplet state 
by ~ 25 meV/T while leaving the energy of the singlet state 
unaltered so that it can be readily accounted for later, if 
necessary. 

The diagonalization procedure for the two-particle Hamil- 
tonian is performed by expanding two-electron spinless wave 
function in the basis 

*s(ri,r a ) = £/3tf hMri)i&;(r a ) + (-l) S ^(riM(r 2 )], 

ij 

(4) 

which is symmetric for singlet (S = 0) state and antisymmetric 
for triplet (S = 1) state. The summation is carried over i < j 
for the singlet and i < j for the triplet. 

From the two-particle wave function, we compute the elec- 
tron density as 



Oi) = J |*s(>i,r 2 )| 2 dr 2 . 



(5) 



The chemical potential of the structure is related to the total 
energy of the system as [6] 

H{Ni + N 2 ) = Eq(N\ + N 2 ) - Eq(N\ +N 2 - 1), (6) 

where E G (N) [note that E G (0) = 0] is the ground state 
energy of the TV-electron state. The exchange energy is given 
by 

J = £g(2)-£g(2), C7) 

where Eq(2) and Eq(2) denote the ground state energy for 
the triplet and singlet state, respectively. 

III. Results 

In Fig. Qa), main panel, we plot the exchange energy as a 
function of the magnetic field for dots with different aspect 
ratio r = R y /R x , R x = 30 nm. In insets (I), (II), and 
(III), the contour plots of the two-dimensional potential are 
displayed for r = 1, r = 1.5, and r = 0.67, respectively. 
Other parameters for the model potential are Vl = Vr = 25 
meV, d = 60 nm. It is observed that the magnitudes of 
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Fig. 1. (a) Main panel: exchange energy dependence on the magnetic field 
for dots with different aspect ratios r = R y /R x , Rx = 30 nm. The insets 
(I), (II) and (III) show the contour plots of the confining potential for r = 1, 
r = 1.5 and r = 0.67, respectively, (b) Exchange energy maximum [J(B=0)] 
and magnetic field at singlet- triplet transition [B(J=0)] and (c) merit factor 
J/hu) c [w c = eB/m t c) as a function of the aspect ratio r = R y /R x - For 
(a), (b) and (c), Vl = Vr = 25 meV, d = 60 nm. 



maximum and minimum values of the exchange energy are 
larger for a larger aspect ratio r because of enhanced inter-dot 
coupling. As the aspect ratio r increases, the magnetic fields 
at the singlet-triplet transition point, at the exchange energy 
minimum, and at the saturation point of exchange energy 
all shift to smaller values, leading to a more compressed 
appearance of the curve in the horizontal direction. In Fig. 
Qb), as the aspect ratio increases from 0.67 to 1.5, the 
exchange energy maximum [ J(B = 0)] ramps up slowly from 
9.5 to 13 /ieV, while the magnetic field at the singlet- triplet 
transition point [B( J = 0)] decreases more rapidly from 4 to 
1.75 T. For spin-based quantum computing, a larger J(B = 0) 
and a smaller B(J = 0) are both desirable because they offer 
better control of the quantum logic gate [1], [8]. 

Therefore, we define a merit factor, J(B = 0)/T%uj c , where 
lo c = eB/m*c is the cyclotron frequency. In Fig. [flc), we 
observe that the merit factor increases almost linearly with 
the aspect ratio, with an enhancement factor of more than 
twice, as the aspect ratio increases from 0.67 to 1.5. It is also 
interesting to note in Fig. ^a) that at B w 0.8 T, the values 
of the exchange energy for dots of different aspect ratios are 
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Fig. 2. (a) Schematic of the charge stability diagram of coupled quantum 
dot systems illustrating the charging process of the first and second electrons, 
(b) Contour plots of chemical potentials (lower branches) and /i s (2) 

(upper branches) at r = R y /Rx = 0.67 (dotted), 1 (dashed) and 1.5 (solid). 
d = 60 nm, R x = 30 nm, B = T, reference value fJ.(l) = £t(2) = —17 
meV. 



very close to each other (J w 8.5 /xeV), which implies that 
at this magnetic field geometric factor plays a minor role in 
determining the exchange energy. 

Figure |2ja) shows the schematic of the charge stability 
diagram of coupled quantum dot systems containing up to 
two electrons [6]. Vl (Vr) denotes the controlling gate bias 
for the left (right) dot. The lower curved branch (solid line) 
corresponds to the bias condition under which the chemical 
potential of the first electron in the quantum dot [/^(l), shown 
by the solid line inside the confining potential] is aligned with 
the source and drain chemical potentials /is and /i£> and the 
first electron enters the system. Here, we consider the situation 
close to equilibrium, so in a first approximation we assume 
Ms — Md- When Vl and Vr are increased, the chemical 
potential for the second electron [/i(2), shown by the dotted 
line inside the confining potential] are lowered and aligned 



with fig = fijj and the second electron is charged into the 
coupled dots. The upper curved branch (solid line) corresponds 
to this situation. The separation along the main diagonal 
between the lower and upper branches is called the double- 
triple point (DTP) separation (also called the anti-crossing 
separation) [6], [7], [12]. The number in the parentheses on 
the left (right) gives the electron number on the left (right) dot 
in distinct stable charge regions. 

In our simulation, we choose a fixed reference value of the 
source and drain chemical potentials = /i S (2) = —17 

meV], and plot on the same diagram the chemical potential 
contours for the first and second electrons with respect to 
two controlling gate biases so as to obtain the charge stability 
diagram. Figure[2jb) shows the chemical potential contours for 
the first electron /i(l) (lower branches) and for two electron 
singlet /i S (2) (upper branches) for d — 60 nm and r = 
Ry/R x = 0.67 (dotted), 1 (dashed), and 1.5 (solid) with 
R x fixed at 30 nm at zero magnetic field. As the aspect 
ratio decreases, the single-particle eigenenergies increase as a 
result of stronger confinement in each dot, hence the chemical 
potential contours shift from the lower left corner to the upper 
right corner. The DTP separation is AVl = AVr = 2.83, 
2.94, and 3.04 meV for r = 0.67, 1, and 1.5, respectively. 
According to the "classical" theory in Ref. [6], a larger DTP 
separation signifies a stronger inter-dot coupling strength. 
Therefore, the inter-dot coupling strength is larger for a larger 
aspect ratio, which is consistent with the above result of a 
larger exchange energy value at zero magnetic field found for 
r > 1 [see Fig. [Tla)]. 

In Fig. |3] the exchange energy is shown as a function of Vl 
at Vr — 21 meV, or in other words, as a function of inter-dot 
detuning e = Vl — Vr for coupled circular dots (R x = R y = 
30 nm, d = 60 nm) at different magnetic fields. The exchange 
energy at different magnetic fields is almost invariant over the 
V L interval from 21 to 29 meV. But for V L > 29 meV, it 
increases superlinearly for magnetic field < B < 5 T and 
decreases superlinearly at B = 6 T. For B — 5 and 6 T, the 
J saturation is seen as well. The onset of superlinearity is 
consistent with the difference in the electron density between 
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Fig. 4. Exchange energy J as a function of Vl in the model potential at 
different magnetic fields. Vr is fixed at 21 meV. Other parameters are d = 60 
nm, R x = 30 nm, R y = 45 nm. For clarity, curves for different B fields are 
shifted vertically by —0.3 meV. 



the singlet and triplet states. As shown below in column (II) 
of Fig. |5J as Vl becomes larger, the electron density of the 
singlet state becomes localized in the left dot while the electron 
density of the triplet state continues to be spread over the two 
dots. The difference in electron density and the corresponding 
elecron-electron interaction becomes larger with the increasing 
inter-dot detuning, and hence, the magnitude of the exchange 
energy increases. 

Figure |4] shows the dependence of the exchange energy on 
Vl at Vr = 21 meV, for coupled elliptic dots (R x = 30 
nm, R y = 45 nm, g? = 60 nm) at different magnetic fields. 
The qualitative behavior of the exchange energy for r = 1.5 is 
similiar to that for r = 1 in Fig.|5J except that saturation occurs 
for all the investigated magnetic fields at smaller inter-dot 
detuning. In general, the onset of exchange energy saturation 
occurs when the inter-dot detuning is so large that the electrons 
in both singlet and triplet states starts to localize in the same 
dot and the difference in the electron density between these 
states becomes smaller. We observe that exchange energy 
saturation occurs for Vl > 35 (e = Vl — Vr > 14) meV for 
r — 1 in the absence of the magnetic field (not shown). This 
detuning value for the onset of exchange energy saturation is 
larger than that for the r = 1.5 case because for r = 1.5 the 
effective size of each dot is larger and the electron density 
for both the singlet and triplet states more easily localize into 
a single dot. The general trend in the superlinear increase of 
the exchange energy as a function of the inter-dot detuning 
is consistent with recent experimental observations [8]. For 
d = 60 nm, the exchange energy in the saturation regime is 
in the meV range, which is typical of a single dot of the same 
size (R = 30 nm) [4], [11] and two orders of magnitude larger 
than the value at e = 0. 

In Fig. |3 we show the two-dimensional electron density 
plots for (a) r = R y /R x = 1 and (b) r = R y /R x = 1.5. In 
each case rows (I) and (II) are for the singlet and triplet states, 
respectively. Different columns are for different parameters: (I) 
V L = V R ; = 21 meV, B = 0, (II) V L = 29, V R = 21 meV, 
B = and (III) V L = 29, V R = 21 meV, B = 6 T. In all cases, 
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Fig. 5. Two-dimensional electron density plots for (a) r = R y /R x = 1 
and (b) r = R y /R x = 1.5, In each case rows (I) and (II) are for the singlet 
and triplet states, respectively. Different columns are for different parameters: 
(I) V L = V R = 21 meV, B = 0, (II) Vl = 29, V R = 21 meV, B = 
and (III) V L = 29, V R = 21 meV, B = 6 T. For all cases, d = 60 nm and 
Rx = 30 nm. 



d = 60 nm and R x = 30 nm. By comparing the data in Fig. [5] 
(a) and (b), we see that at zero interdot detuning (Vl — Vr) 
and B = [column (I)], the electron density of both the singlet 
and triplet states shows two peaks localized at the center of 
each dot which are symmetric with respect to x = 0. The 
exchange energies for in these cases are comparable to each 
other (J = 15 and J = 19 fisV, respectively). At non-zero 
interdot detuning (Vl — Vr — 8 meV) and B = [column (II)], 
the singlet electron density in (a) for r = 1 shows a secondary 
peak in the right dot, while this peak is absent for a larger 
aspect ratio r = 1.5 in (b). This is because for a larger aspect 
ratio, as was already mentioned above, the onset value of the 
inter-dot detuning for two electron localization into a single 
dot is smaller (see discussion on Figs. [5] and |4] above) due to 
the relaxation of the confinement potential in the y direction. 
In this case, the triplet electron density profiles for both cases 
show a secondary peak in the right dot, and are similar to one 
another while the calculated exchange energies are J = 297 
and J = 1050 /ieV, respectively. At non-zero interdot detuning 
(V L -Vr = 8 meV) and B = 6 T [column (III)], there is a 
drastic difference in the electron density profiles between the 
results in (a) and (b): in (a) for both singlet and triplet the 
electron density shows a higher (lower) peak localized at the 
center of the left (right) dot, while in (b), for both singlet and 
triplet the electron density is totally localized in the left dot 
with two equal peaks at y = ±8 nm, x — —30 nm. Due to the 
large overlap of the electron wavefunctions in the latter case, 
the exchange interaction J = —577 /ieV is more than two 
orders of magnitude larger than J = — 1 /ieV in the former 
case. 
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IV. Conclusion 

We have computed the exchange energy and derived the 
stability diagram for model coupled elliptic quantum dot 
systems with different aspect ratios. We find larger exchange 
energy and smaller magnetic field at singlet-triple transition for 
dots with larger aspect ratios R y /R x > 1, which provides a 
better control of the exchange interaction. By investigating the 
double-triple point separation and the curvature of chemical 
potential contour lines on the stability diagram for different dot 
configurations, we find the inter-dot coupling decreases with 
increasing inter-dot separations, increasing magnetic fields, 
and decreasing dot aspect ratios. In the weak coupling regime 
(inter-dot distance 60 nm), we find three distinctive regions in 
the dependence of the exchange energy (J) on the inter-dot 
detuning (e): (1) For small e, J maintains a nearly constant 
value and the electron density of the singlet and triplet states 
are spread between the two dots and their difference is small. 
(2) For intermediate e, J depends superlinearly on e. Singlet 
electron density is well localized into one dot, while triplet 
electron density is still spread between the two dots. (3) For 
large e, the increase of J tends to saturate and both singlet 
and triplet electron densities begin to localize into one dot. At 
a fixed aspect ratio, the onset of the superlinear and saturation 
regions depends on the magnetic field. For a constant magnetic 
field, the onset of these two regions occurs at a smaller inter- 
dot detuning value for dots with larger aspect ratios. 
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